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A computation in adiabatic quantum computing is implemented by traversing a path ol nonde- 
generate eigenstates ol a continuous lamily ol Hamiltonians. We introduce a method that traverses 
a discretized form ol the path: At each step we apply the instantaneous Hamiltonian for a random 
time. The resulting decoherence approximates a projective measurement onto the desired eigenstate, 
achieving a version ol the quantum Zeno effect. II negative evolution times can be implemented with 
constant overhead, then the average absolute evolution time required by our method is 0(L 2 /A) for 
constant error probability, where L is the length ol the path ol eigenstates and A is the minimum 
spectral gap ol the Hamiltonian. The dependence of the cost on A is optimal. Making explicit the 
dependence on the path length is useful for cases where L is much less than the general bound. 
The complexity of our method has a logarithmic improvement over previous algorithms of this type. 
The same cost applies to the discrete-time case, where a family of unitary operators is given and 
each unitary and its inverse can be used. Restriction to positive evolution times incurs an error 
that decreases exponentially with the cost. Applications of this method to unstructured search 
and quantum sampling are considered. In particular, we discuss the quantum simulated annealing 
algorithm for solving combinatorial optimization problems. This algorithm provides a quadratic 
speed-up in the gap ol the stochastic matrix over its classical counterpart implemented via Markov 
chain Monte Carlo. 



I. INTRODUCTION AND SUMMARY OF 
RESULTS 

Quantum algorithms are often described by means 
of quantum circuits: the algorithm starts with a well- 
characterized pure state; a sequence of elementary (uni- 
tary) gates is applied; and a final projective measure- 
ment in a fixed basis extracts the result. The circuit 
model may not be best for describing all quantum infor- 
mation processing systems. Adiabatic quantum comput- 
ing (AQC) pQ, sometimes also called quantum anneal- 
ing [2j |3j SJ [5j [6] , has been proposed as an alternative. 

In AQC the computation is performed by smoothly 
changing the interaction parameters of the Hamiltonian 
under which the system evolves. The initial state is a 
nondegenerate eigenstate of the Hamiltonian. The adia- 
batic theorem of quantum mechanics asserts that if the 
continuously related eigenstates remain nondegenerate 
and the Hamiltonians change sufficiently slowly, then the 
final state of the system is close to the continuously re- 
lated eigenstate of the final Hamiltonian [7] . The last step 
is a standard projective measurement. AQC is polyno- 
mially equivalent to the quantum circuit model 8 . 

In this paper we give a method for traversing eigen- 
state paths of Hamiltonians that differs from AQC by the 
use of evolution randomization. The method is based on 
previous results [9j [10] in which the evolution of AQC is 
replaced by a sequence of projective measurements onto 
the instantaneous eigenstate of the Hamiltonian with the 
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phase estimation algorithm |llj , which exploits the quan- 
tum Zeno effect. Both AQC and the Zeno-based model 
work, in essence, because an effective level decoupling is 
introduced in the Hamiltonian eigenbasis by phase can- 
cellation in AQC or projections in the Zeno case. Our 
method also implements a version of the quantum Zeno 
effect. We choose a discretization of the eigenstate path 
and apply the Hamiltonian corresponding to each point 
for a random time. The probability distribution over 
time may be discrete or continuous. Consequently, the 
randomization method can also be used in the case where 
we are given a path of efficiently implementable unitary 
operators and an eigenstate of the last operator on the 
path is to be prepared. This case occurs in the quan- 
tum simulated annealing (QSA) algorithm constructed 
in Ref. [T2]. The probability distribution over evolution 
times must be chosen so as to cancel unwanted coherences 
and simulate the Zeno effect. 

The algorithmic complexity of the randomization 
method is defined as the average sum of the absolute 
evolution times for the Hamiltonians or by the average 
number of times the unitaries are applied. The complex- 
ity can be bounded in terms of a lower bound A on the 
absolute value of the minimum spectral gap of the Hamil- 
tonians or the minimum phase gap of the unitaries, the 
length L of the path of the states (defined below), and 
the desired maximum error e of the final state compared 
to the target eigenstate. We show that the complexity 
is 0(log(L/e) a L 2 /(eA)), where a = if we can evolve 
backward and forward in time, and 1 otherwise. Back- 
ward evolution is possible at the same cost by reversing 
quantum circuits for the forward evolution, if such evolu- 
tion is circuit-based. To achieve this complexity without 
additional dependencies, we use a parametrization of the 
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operators along the path for which the eigenstates move 
at a rate that is close to uniform (up to a constant factor) . 
In many cases of interest, L £ O(l) so that the complex- 
ity is of order l/(eA) up to logarithmic factors. The scal- 
ing with the gap is optimal and is better than the 1/A 3 
of rigorous proofs of the adiabatic theorem [T3j [14j [15] . 

An advantage of our approach is that the only require- 
ment on the Hamiltonians or the unitaries along the path 
is that the length of the desired eigenstate path is well- 
defined. A sufficient condition is that the time derivative 
of the operators exists. In terms of bounds on the Hamil- 
tonians and their derivatives, the worst-case bound is of 
order ||ij|| 2 /(eA 3 ) up to logarithmic factors. This bound 
comes from the inequality L < ||£f||/A [see Eq. (13)] 



and does not depend on a reparametrization. It is bet- 
ter than the known worst-case bounds associated with 
the adiabatic theorem |13j in that it does not depend 
on existence of, or bounds on the second derivative of 
H. The scaling of the bound with the error is worse, in 
that it can be made logarithmic for analytic Hamiltoni- 
ans paths at the cost of less favorable dependencies on 
the other parameters of the problem [TUEI]. Logarith- 
mic scaling can also be obtained with the randomization 
method provided the final eigenstate's energy is known. 
To achieve this scaling, one can use high precision phase 
estimation to determine whether the desired eigenstate 
has been obtained and repeat the algorithm if not. 

Our method is intuitively explained by the quan- 
tum Zeno effect. Suppose that the path of states 
\ip(li)), . . . , \ip(l q )) satisfies the condition that for each 
j, \4>(lj-i)) is sufficiently close to If we initialize 

the state \ip(h)} an d sequentially apply projections onto 
the we prepare \ip(l q )) with good probability. We 

first consider an idealized strategy, where the projections 
are replaced by quantum operations of the form 

M h (p) = P ljP P h + £((1 - P h )p(l - P h )) , (1) 

where Pi = and £ is an arbitrary quan- 

tum operation that may vary from instance to instance. 
This can be thought of as a projective measurement of 
p onto \ip(lj)) followed by a process that does not af- 
fect The fundamental effect of Mi. is to remove 
coherences between \ip(lj)) and orthogonal states. It is 
this decoherence that induces the quantum Zeno effect by 
suppressing transfer of population to orthogonal states. 
An approximation of this effect is achieved if we replace 
the Mi . by random applications of Hamiltonians or uni- 
taries with \ip(lj)) as an eigenstate. We formalize this 
claim in Sec. |TT] Thm. [T] and give an upper bound for 
the error in the approximation in terms of the character- 
istic function of the probability distribution underlying 
the randomization. 

We focus on the Hamiltonian-based version of the ran- 
domization method. The analysis for the unitary version 
is a straightforward discretization. In the Hamiltonian 
version, the randomization method takes as input a con- 
tinuous path of Hamiltonians H — {H(s), s E [0, 1]}, and 



a nondegenerate eigenstate \ip(0)) of H(0). The method 
aims to output the corresponding nondegenerate eigen- 
state of -ff(l), denoted by \tp(l)}, with high fidelity. 

We require that the eigenstates |^(s)) are nondegen- 
erate with A a lower bound on the energy gap. If \tp(s)} 
is differentiable (see Appendix VI for the more general 



case), we can assume without loss of generality that the 
phases of the \ip(s)) are chosen geometrically, so that 
(d s ip(s) \i>(s)) = 0, which gives a path length 



L 



\dsi>(s)) \\ds . 



(2) 



The quadratic cost dependence on L comes from a simple 
Zeno effect when an ideal decoherence process accord- 
ing to Eq. ([I]) is used. It is probably not fundamental: 
Coherent versions of the adiabatic path achieve scalings 
O(L) [12]. The dependence of the cost on 1 /A is unavoid- 
able for methods with only oracle access to the Hamilto- 
nian or unitaries. This can be seen intuitively by noting 
that we must, in a sense, distinguish between the desired 
eigenstate and the others, which requires that we evolve 
the relative phases sufficiently far. More rigorously, an 
asymptotically better dependence would result in an un- 
structured search algorithm better than Grover's, which 



is known to be impossible. See Sec. IV A 



The paper is organized as follows. In Sec. [H] we ex- 
plain how the quantum Zeno effect can be exploited, show 
how to approximate projective measurement operations 
by means of evolution randomization, and discuss sev- 
eral probability distributions that are useful for random- 
ization. The randomization method and its complexity 
are analyzed in Sec. |III[ In Sec. |IV A| we show that the 
randomization method provides the expected quadratic 
quantum speed-up for the unstructured search problem. 



In Sec. IV B we describe the QSA to simulate slowly vary- 



ing classical Markov chains. In Sec.[V]we show the equiv- 
alence of our randomization method with a coherent ver- 
sion of the quantum Zeno method implemented via the 
phase estimation algorithm, and briefly discuss related 
works. We summarize in Sec. IVII 



II. RANDOMIZED EVOLUTIONS 



A. Adiabatic quantum computing using the Zeno 
effect 



The quantum Zeno effect is based on the fact that, 
for a small displacement 5', the probability of project- 
ing \tp(s + 5')) onto \ip(s)) decreases with (<5') 2 , while the 
distance between states is linear in 5' [T8j [19l 12Q] ■ There- 
fore, for the path of states {\i{j(s)}}, the final state IV'(l)) 
can be prepared from the initial state IV'(O)) with high 
fidelity by use of a sequence of measurement projections 
onto intermediate states \ip(si)},--- ,|^(s g )), < si < 
• ■ ■ < s q = 1 . We choose Sj so that the fidelity of the 
final state with respect to | -0 ( 1 ) ) is sufficiently close to 
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unity. It is not necessary to keep track of the measure- 
ment results at intermediate steps, which gives rise to the 
following definition. 

Definition 1. A projective-measurement operation onto 
\ip{l)) is a quantum operation of the form 

M,(p) = P lP Pt + £((t - P0p(l - P t )) , 

with Pi = \ip(l)}{i^(l)\ and £ arbitrary quantum opera- 
tions that may vary with I. 

We assume a monotonically increasing parametrization 
s(l), with I e [Q,L'], s(0) = and s(L') = 1. We define 

1^(0) = IVK s (0))- (Objects with a tilde correspond to 
objects in the new parametrization). Later we consider 
s(l) so that II = L, the path length of Eq. We for- 
mulate the Zeno method for quantum state preparation 
as follows [9j[10j EE] : 

Lemma 1 (Zeno effect). Consider a continuous path of 
states {\ip(l))}ie.[o L'] and assume that, for fixed d and all 
S, 

\(mm+6))\ 2 >i-d 2 s 2 . 

Then the state \^(L')) can be prepared from \ip(fi)) with 
fidelity p > by \{L') 2 d 2 /{I— p)] intermediate projective- 
measurement operations. 

Proof. Divide [0, L'] into q = \(L') 2 d 2 /(l -p)] equal seg- 
ments and set S = L' jq. At every point lj = jS, 1 < j < 
q, we perform a projective-measurement operation onto 
\^{lj)). The final state is M iq o Mi q _ 1 o • • • o M h (p), with 
p = \ip(0)) (ip(0)\. The output fidelity is bounded as 

tr [P lq (M lq o • • • o M h (p))] >\\P q --- P h \m) || 2 

= ^ =1 m ] )\^ 1 ))\ 2 

> (1 - d 2 5 2 ) q > 1 - d 2 L' 2 /q >l-(l-p)=p (3) 

From Lemma [T] and assuming a uniform parametriza- 
tion, defined to satisfy L(s(l)) = L(l) = I, d = 1, and 
L' = L (see Appendix |VI| , it follows that the state IV'(l)) 
can be obtained with fidelity p starting from | -0(0)) with 
0(L 2 j(\ — p)) projective-measurement operations. 

B. Approximating projective-measurement 
operations through randomized evolutions 

We assume that evolutions under H(s) for time t can 
be implemented at a cost linear in |i|||i/(s)||, as in AQC. 
That is, we do not take into account the cost of simu- 
lating H(s) for small time intervals. By rescaling H(s) 
if necessary, we can assume that ||iJ(s)|| < 1. Thus, the 
cost of the randomization method is determined by the 
sum of the absolute evolution times. Although we con- 
sider the case where the evolution time t can be negative, 



one often restricts t to be nonnegative. This restriction 
is justified if the Hamiltonians are physical without a 
simple time-reversal procedure, rather than induced by 
quantum circuits. In the latter case, evolving for nega- 
tive t is as efficient as for positive t and can be realized 
by reversing the quantum circuits. 

We denote by A(s) the spectral gap for the eigenstate 
\ip(s)) of Hamiltonian H(s). The following results also 
apply to the unitary case where we are given operators 
U(s) and A(s) is the phase gap. In the unitary case the 
distributions over time that are used for randomization 
must be concentrated at the integers, and correspond to 
the number of times the unitaries are applied. 

According to Lemma [T] the Zeno method does not re- 
quire that we keep track of intermediate measurement 
results. Thus, any purely dephasing mechanism in the 
instantaneous eigenbasis of H(l) implements a version of 
Mi. A natural choice for such a decoherence mechanism 
is the evolution induced by H(l) for a (unknown) ran- 
dom time. This is the subject of next theorem, where we 
bound the residual coherences in terms of the character- 
istic function of the random time distribution. 

Theorem 1 (Randomized dephasing). Let \ip(l)) be a 
nondegenerate eigenstate of H{l), and {u>j} be the en- 
ergy differences to the other eigenstates \ipj(l))- Let T be 
a random variable associated with the time of evolution 
under H(l), and IZf the corresponding quantum opera- 
tion. Then there exists a quantum operation £ such that, 
for all states p, 

\\{Mf -Kj){p)\\ tI <e = sup 1*^)1 , 

UJj 

where Mf is the projective-measurement operation de- 
fined in Definition^ with £ specified, and $ is the char- 
acteristic function of T. 

We give the proof in Appendix |VI| It is based on com- 
puting the coherences after the randomized evolution in 
terms of the characteristic function of T as 

nf(\m)(^M) = H^)\m)(Mi)\- (4) 

The average cost of randomization is given in terms of 
the random variable T as (|T|), the expected value of the 
absolute evolution time. If T takes only positive values 
the average cost is given by 

(T) = -i[d u $}(0) , 

provided it is finite. Note that if T < is allowed, then 
the average cost can be reduced by shifting T's distribu- 
tion so that is a median of T. 

We can bound the required average cost per step from 
below by f2(l/A), with A a lower bound on the smallest 
gap inf s |A(s)|, by means of the following theorem: 

Theorem 2. Let T be a random variable with charac- 
teristic function Then, for all u>, 

cost(T) = (|T|) > . 
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Proof. From the definition of we obtain 

l-|$H|<|l-$(o;)|< J |1 - e*"*|d/i(i) 

< / \ut\diJt) = (\T\)\u\ , (5) 



with /i the probability distribution of T 

We want to ensure that after the randomized evolution, 
the remaining coherences bounded by |$(w)| for \u>\ > 
A are small. Because of Thm. [2] the average absolute 
evolution time can be bounded by f2(l/A). 

If T < is permitted, the bound of Thm. [2] can be 
achieved up to a constant factor. See Example [2] and 
Lemma [3] in Sec. |II C| For the case where we are given 
a path of unitaries and T is restricted to the integers, 
it suffices to consider the characteristic function on the 
interval [— n, tt]. The results of this section are otherwise 
unchanged. 

Repetition of the randomized evolution step decreases 
the error exponentially fast in the number of repetitions, 
as shown by the following argument. For independent 
random variables T± and T 2 , the characteristic function 
of the sum T' = T x + T 2 is 



(6) 



with <&i the characteristic function of X^. Thus we have 
the following lemma (the notation is that of Thm. [IJ: 

Lemma 2. Let T be a random variable with character- 
istic function <I>, and sup Wj |<£(wj)| = e. Let T" be the 
sum of n independent instances of T. Then there exists 
a quantum operation £ such that for all states p, 



\(M( 



){p)\\ 



< e" 



C. Examples of randomized evolutions 

We consider some examples of randomized evolution 
steps involving different time distributions. 

1. Consider the case where all the orthogonal eigen- 
states to \ip(l)) are degenerate and the spectral gap 
of H(l), denoted by ui±, is known. We then choose 
a random variable T UJl that takes the values t = 
or t = tt/lui, each with probability 1/2. The aver- 
age cost is it/(2uji). The characteristic function for 
this distribution satisfies 



|*MI 



cos 



TTUJ 
2t^L 



(7) 



Since 3?(oji) = , Thm. [T^implies that the pro- 
jective measurement onto \ip(l)) can be simulated 
exactly with this distribution. The assumptions in 
this example may seem unrealistic, but it provides a 
basis for the randomization method in unstructured 



method to the case where the spectrum is known. If 
there are k distinct absolute eigenvalue differences 
ujj, the independent sum of T u . has the property 
that the characteristic function is identically zero 
on the eigenvalue differences. The average cost is 

2. Let T's probability density be proportional to 
sinc(At) 4 , A > 0. The function sine is defined 
as sinc(i) = sin(i)/i. The Fourier transform of 
Asinc(Ai) /it is the indicator function of the interval 
[—A, A]. The characteristic function of T is there- 
fore proportional to the four-fold convolution of this 
indicator function with itself, which is continuous 
and has support [— 4A,4A]. There is no error in ap- 
proximating the projective-measurement operation 
by randomized evolution if we choose 4 A = A, with 
A a lower bound on the minimum gap. The aver- 
age cost (\T\) is proportional to 1/A = C(l/A). 
According to Thm. [2] this is optimal. A possible 
problem is that the tail distribution of T is large: 
Moments of order greater than 2 are unbounded. 
Lemma [3] shows that this can be remedied. For the 
unitary case we modify T by restricting to the inte- 
gers. That is, we set Prob(T = n) cx sinc(nA) 4 . For 
A < 7r/4, the restriction of the characteristic func- 
tion to [— 7r,7r] is unchanged (see Lemma [1]), so for 
the case where the cigenstate path is determined 
by a path of unitary operators, the same average 
cost of 0(1/ A) is obtained. 

3. When the eigenstate path is determined by unitary 
operators, a simple choice of T is the uniform dis- 
tribution on integers between and Q — 1, where 
Q = [2-7r/A~|. If we repeat the randomization step 
n times, we can bound the error with respect to the 
desired projection by (Lemma |2| 



SUp \$(Uj)\ 



< 



1 1 



Q 1 



oi'A 



< 



1 

2" 



(8) 



The average cost is n(Q- l)/2 G 0(n/A). To have 
error at most e, the cost is C(log(l/e)/A). 

If negative T can be used, we can shift T by 
— LQ/2J . This does not affect the absolute values of 
the characteristic function but reduces the average 
cost by a constant factor near 1/2. 

4. If T is unrestricted, we can consider T with Gaus- 
sian distribution Af(0, a). Note that restricting to 
0-mean Gaussians minimizes (ITI) since the mean 
and the median coincide. The absolute value of the 
characteristic function is 



(9) 



=exp - 



search (Sec. IV A I. It is possible to generalize the 



The error of the randomization step with respect 
to the desired projection is bounded by 

,2 A 2 , 



|*(A)| = exp (-(T 2 A 2 /2) 



(10) 
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For this distribution, (\T\) = o-y/2/w. To have er- 
ror at most e, we need a > 21og(l/e) 1//2 /A. This 
gives an average cost of C[log(l/e) 1 / 2 /A]. 

If T must be positive, we can displace the Gaussian 
by x > and condition on positive outcomes. The 
error can be estimated as the sum of the probability 
that the Gaussian is negative, which is bounded by 



-*7(20 



, and the right-hand-side of Eq. ( 10 1 . The 



average cost is 0(x + a). To have error at most e, 
let a = 21og(2/e) 1 /2/A an d x = ^er log(2/e) 1 / 2 . 
The average cost is then C[log(l/e)/A]. According 
to Thm. [3] (below) this is optimal for positive T. 

5. In the case where T must be supported on inte- 
gers, one can try to approximate the Gaussian by 
the shifted binomial distribution obtained from the 
sum of 2m independent { — 1/2,1/2} mean-0 ran- 
dom variables. The absolute value of the charac- 
teristic function is 



|<&MIHcos(u,/2)| s 



(11) 



This requires m € 6(log(l/e)/A 2 )) to achieve 
error e in approximating the desired projection. 
The average cost is then C(log(l/e) 1/2 /A)). As 
in Example HJ we can shift the distribution by 
0[log(l/e) 1 / 2 /A] and condition on positive integers 
to ensure that T is positive and obtain an average 
cost of C[log(l/e)/A]. 

Except for Example [2j the distributions above do not 
achieve the optimal asymptotic cost for unconstrained 
T. In the case of Example [2j the probability density 
determined by sinc(Af) 4 has long tails and unbounded 
moments. This is improved by the following lemma. 

Lemma 3. There exist probability densities f for T that 
achieve cost (\T\) = 0(1/A) and error |$(w)| = for 
\uj\ > A, and that have bounded moments of all order, 
i.e. ((T — (T)) n ) < 00 V n > 0. 

We give a constructive proof in Appendix |VI| For these 
distributions and A < tt, discretization does not result 
in an increase in the error, see the next lemma. Note 
that in the discretized case we are only interested in the 
region of eigenphases [—tt, tt] and the relevant gap is the 
eigenphase gap. 

Lemma 4. Let f be a probability density whose char- 
acteristic function has support in (—A, A) with A < tt. 
Then the restriction of f to the integers is a well-defined 
probability distribution with prob(k) — f(k) and charac- 
teristic function <t>(c<j) = for \oj\ £ (A,7r]. 



We give the proof in Appendix |VI) 

For positive T and if only a lower bound A on the gap 
is known, it is not possible to improve asymptotically 
over the shifted and conditioned Gaussian distribution 
of Example |4j 



Theorem 3. Let T be a positive random variable with 
characteristic function Then sup| w | >A |$(w)| > 

The proof is in Appendix |VI| 



III. THE RANDOMIZATION METHOD 

The goal of the randomization method is to prepare 
the nondcgcnerate eigenstate IV'(l)) °f -^(1) by travers- 
ing of the path \ip(s)}. This path is determined by 
the family H = {H(s)}. Ideally, we choose the uni- 
form parametrization s(l) discussed in Sec. |II A and 
Appendix |VI| Under such a parametrization the cigcn- 
states \ip(s(l))) move at a constant unit rate along the 
path. Finding the uniform parametrization is difficult 
in general. We therefore consider an arbitrary subuni- 
form parametrization I 6 [0, L'] 1— > s(l) so that the rate 
at which the states move is bounded by unity. Note 
that L' > L, with L the path length. A subuniform 
parametrization can usually be obtained from known 
properties of H(s); see Lemma [5] and Eq. (14) below. 
We discretize the path using q £ 0((L') 2 ) segments in 
order to achieve bounded error. 

The randomization method uses randomized evolu- 
tions TZf to approximate the projective-measurement op- 
erations Mi at values s(l). Here, I = kd for k = 1, . . . , q, 
with q — L' /S and S sufficiently small. For good asymp- 
totic behavior, we choose T as in Lemma|3]or Example|2] 
If T must be positive, we use the shifted and conditioned 
Gaussian distribution of Example |4j If T must be re- 
stricted to the integers, as in the case of a path U(s) of 
unitaries, we use the discretized version of T (Lemma EJ. 
We obtain: 

Theorem 4 (Randomization method). There are 
choices of q and T in the randomization method such 
that the method outputs |^(1)) starting from 1*0(0)) with 
fidelity at least p and average cost 



O 



(L'f(\og(L'/{\-p))y 



(1-P)A 

where a = if T can be negative and a — 1 otherwise. 

Proof. We choose a step increment S — L'/q, with 
q = [2(L') 2 /(1 — p)~|. For this choice, Lemma [l] guar- 
antees that, if we were to implement the projective- 
measurement operations exactly, the error in the prepa- 
ration of IV'(l)) would be bounded by (1 —p)/2, because 
d < 1 for subuniform paramctrizations. Wc need to 
choose T such that the additional contribution to the 
error due to the differences between the randomized evo- 
lutions and the projective-measurement operations is also 
bounded by (1 —p)/2. Suppose that the error according 
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to Thm. [T] is bounded by e. After r steps we have 

\\(M lr o...oM h -Rlo...oRf i )(p)\\ ii 
= \\(M lr o • • • o M h - R lr o M lr l o • • • o M h ){p) 
+ {Rj r o M lr _ t o...oM h -Rlo...o Rl){p)\ 

< \\(M lr -R lr )(a)\\ tI 
Ri r ( M i r -i o-.-oM h - Rf ri o • • • o Rl)(p) 

< e + (r — l)e = re , 



tr 

(12) 



where we used the fact that quantum operations are 
trace-norm contracting, and we implicitly applied induc- 
tion in the last steps. The desired bound on the error 
requires e < (1 -p)/(2q) = (1 - p) 2 /(4(L') 2 ). According 
to Lemma [3] and Example [4] this can be achieved at an 
average cost (|T|) of 0(1/A) if T can be negative, and 
C(log(l/e)/A) otherwise. The total cost for the proce- 
dure is 0[qlog(q/(l — p)) a / A], and substitution of the 
value for q yields the claimed bound 

For differentiable H(s) and eigenstate path \i()(s)), 
we can obtain a subuniform parametrization s(l) from 
bounds on the derivative of H(s) and the gaps. For this 
we need the next lemma. 

Lemma 5. Suppose that H(s) is differentiable and 
HV'(s))} is a path of nondegenerate eigenstates of{H(s)} 
with spectral gap A(s) =/= 0. Then 



\ds*P(s)) || < 



djm ii 

|A( S )| 



The proof is in Appendix |VI| 

Define \\H\\ = sup s \\d s H(s)\\. We obtain 



\d.i/j(a)) \\ds < V 



\H\\ 
~A~ ' 



(13) 



with A a lower bound to the minimum absolute value of 
the gap. This L' is achieved for the parametrization 



s(l) = , 
ll#ll 



(14) 



which is subuniform in general. Using this parametriza- 
tion we obtain the following corollary: 



To conclude this section we consider the following two 
questions: What is the probability that the cost of the 
randomization method exceeds the average cost by a con- 
stant factor? How does the actual path followed by the 
states obtained in a given instance of the randomization 
compare to the adiabatic path? 

The average cost of the randomization method is 
(C) = q(\T\), where q is defined in the proof of Thm. [4] 
with T the relevant random variable. The probability 
prob(C > a(C)) is therefore at most 1/a (Markov's in- 
equality). If the higher-order moments of T are bounded, 
better bounds can be obtained. In particular, for the dis- 
tributions whose characteristic functions have smooth, 
compact support, prob(C > a(C)) decreases superpoly- 
nomially in a. For T based on Gaussians, the decrease 
is e' n< - a2 \ Since C is determined by a sum of q inde- 
pendent instances of |T|, better bounds can be obtained 
for specific choices of T, particularly if q is large. In par- 
ticular the variance of C is inversely proportional to q if 
T has finite variance and Chebyshev's inequality or, for 
sufficiently well-behaved T, large-deviation theory can be 
applied. 

A distinguishing feature of the randomization method 
is that any given instance involves unitary evolution, 
which means that the sequence of states obtained is pure. 
What is the probability (over the randomization of the 
evolution times) that every state in the sequence of pure 
states has fidelity at least 1 — 7 with respect to the cor- 
responding eigenstate along the adiabatic path? In view 
of the proof of Thm. [4] the probability that the state 
after the r'th step has fidelity at least 1 — kre with re- 
spect to \ip(l r )) is at least l/k (by Markov's inequality). 
In particular the fidelity of the last state obtained is at 
least 1 — k(l — p) with respect to 1^(1)) with probabil- 
ity l/k. One can deduce that many of the states ob- 
tained in a typical instance of the randomization method 
are close to the corresponding states along the adiabatic 
path. Given that the deviation from the adiabatic path 
executes a kind of random walk, it is reasonable to con- 
jecture that for appropriate choices of parameters, the 
probability that all states obtained are close to the adia- 
batic path is also high. 



IV. EXAMPLES OF QUANTUM 
COMPUTATIONS VIA EVOLUTION 
RANDOMIZATION 



Corollary 1. Let H(s) be a differentiable path of Hamil- 
tonians and A a lower bound on the minimum absolute 
value of the spectral gap. Then we can prepare IV'(l)) 
from IV'(Q)) with bounded error probability at cost 



° "|F(lo S (ll^ll/A 



where a = if we can evolve for negative times and a = 1 
otherwise. 



A. Unstructured search 

In Grover's algorithm [2Tj we want to find a single 
marked clement S in a space of N — 2™ elements. For 
this, we build the Hamiltonian 

H(s) = -[s\S){S\ + (l- s )|+)(+|], (15) 

acting on a set of n qubits. Here, |+) is the equal su- 
perposition state and \S) the solution state, which is the 
computational basis state corresponding to the marked 
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element. Evolving with H(s) for time t can be done us- 
ing 0(|i| 1+r ') conventional oracle calls, with 77 > ar- 
bitrarily small [22 . For any s, H(s) is nondegenerate 
in the subspace spanned by {|+) , \S}}. If \tp(s)} is the 
eigenstate with largest eigenvalue, we seek to prepare 
IV'(I)) = \S) from \ip(0)) = |+) with sufficiently high 
probability. Preparation of |t/>(1)) using AQC was stud- 
ied in Ref. [53] . 

The energy gap of H(s) can be obtained exactly in the 
relevant subspace. It is 



A(s) = y/l-4s(l- s)(l-l/N) 



(16) 



which is minimized at s = 1/2, giving A = A(l/2) = 
1/s/N. The path length L can also be obtained exactly 
and, for large N, we have L w tt/2 (the states |+) and \S) 
are almost orthogonal) . From Thm. [4] the average cost 
of the randomization method for constant probability of 
success is 0(1/ A) £ 0(l/\f~N) if the parametrization 
is uniform. In the large N limit, this parametrization 
satisfies 



us 1 I 1 1 - s 

Us) ~ - arctan — — — ; 

V ; 2 \s/Nl/2-s 



(17) 



which satisfies < l(s) < tt/2. The randomiza- 
tion method then consists of a sequence of projective- 
measurements operations at values 



1 



cot (21 j) 



2VN 



cot(2j<5) 



2VN 



(18) 



for some S > 0. Note that this is the same evolution path 
as the one considered in Ref. [53], and that the rate of 
change of s as a function of I is A(s(l)). A possible choice 
for S is 7r/4. At s = 1/2, we can implement the phase 
randomization by evolving under if (1/2) for time or 
7r/A, each with probability 1/2. This is the distribution 
in Example I] of Sec. |II C[ and was also used in Ref. [9]. 
It outputs the desired state almost half the time. 

When more than one marked element exist, the above 
randomization method can still be used to output a so- 
lution with bounded error probability: the main effect 
of adding new projectors in H(s) is an increased spec- 
tral gap A'O) > A(s). Thus, the induced decoherence 
still simulates an appropriate measurement in the new 
eigenbasis. If the uniform distribution is used for the 
randomization, then the algorithm is equivalent to the 
one discussed in Ref. [2U, Sec. 8.4. 



B. Quantum simulated annealing 

As the previous example demonstrates, distinguishing 
between the cost induced by the path length and the one 
induced by the gap has important advantages, in par- 
ticular when L E 0(1). Without this distinction, the 
actual cost of the method can be highly overestimated. 
In Ref. |12j we studied quantum simulations of classical 



annealing processes via evolution randomization. An up- 
per bound on the path length in this case is independent 
of the minimum spectral gap T of the classical Markov 
chain (i.e., T is the difference between 1 and the sec- 
ond largest eigenvalue of the stochastic matrix) . Further- 
more, r can be quadratically increased using Szegedy's 
quantum walks |25l 126] . For bounded error probabil- 
ity, the randomization method using these walks has a 
cost 0(1/a/T), where we are disregarding the depen- 
dency on other parameters such as error probability and 
path length. It provides a quantum speed-up with re- 
spect to simulated annealing using Markov Chain Monte 
Carlo methods, where the cost is 0(l/r). Quantum state 
preparation of Gibbs' states using AQC and the Zeno 
method was previously studied in Ref. |10] , but no quan- 
tum speed-up was obtained. Recently, a unitary version 
of the quantum simulated annealing algorithm (QSA), 
that uses Grover's fixed point method, was introduced in 
Ref. [17] . The unitary version improves the dependence 
of the cost of QSA on output fidelity compared to that 
in Ref. [12]. However, the scaling in the gap is the same. 

Basically, QSA is designed to traverse a coherent ver- 
sion of the classical-state path traversed by classical sim- 
ulated annealing. The quantum state path is in a Hilbcrt 
space of dimension corresponding to the size of the classi- 
cal state space. The classical annealing path we consider 
is determined by tt x ((3) = e~^ E ^ / Z((3), where tt x is the 
probability of configuration x in the stationary (Gibbs) 
distribution. E is the associated energy or cost function, 
(3 is the inverse temperature, and Z(0) the partition func- 
tion. The corresponding path in Hilbert space is given by 
the quantum Gibbs states \il>(/3)) — ^2 X \J tt x (/3) \x). Note 
that a measurement in the computational basis samples 
x with probability tt x ((3). Since 

\dp1>(P)) = E« S > - S M)V^/2 \x) , (19) 

X 

we obtain the following lemma. 
Lemma 6. For f3 £ (0, /?/), 

11 \d P m) 11 - *G8)/2 , 

where <j(f3) is the standard deviation of E at inverse tem- 
perature f3. The path length satisfies L < (3fO-/2, with 
a = sup^ a((3). 

If d' is the size of the classical state space and 7 
is the spectral gap of E, then the state \tp(f3f)), for 
(3f = 0((log d')/j), has high probability amplitude in 
the configuration that minimizes E. With this /?/, we 
have 



L e 



a log d' 

7 



(20) 



That L is bounded independently of T is fundamental 
for the success of QSA. Using Szegedy's quantum walks 
we can boost the gap towards 0(y/T) and achieve the 
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desired cost. The details of this procedure are explained 
in Refs. [El [26]. 

The QSA is basically a sequence of steps, each 
constructed to prepare the states \ip(S)}, |^>(2<5)), 
••• , \4>(f3f)} from the initial state \4>(0)), S -C 1. Ac- 
cording to Ref. [15] . these states can be prepared by a 
version of the Zeno effect in which, at each step, the cor- 
responding Szegedy walk is applied a random number of 
times (see Example [3] Sec. II C I. For this distribution the 
cost of the QSA is 

O [ ^logL] e O ( ^^/J ' log(alogd' h) 



/r 



7 



/r 



with r the minimum gap of the Markov chain along the 
path. The results in Sec. |II C| show that using the in- 
verses of the quantum walk steps, the second logarithmic 
factor can be dropped. Because of the way the quantum 
walk is constructed, circuits for the inverses can be ob- 
tained by direct reversal of the circuits for the quantum 
walk steps. In Ref. [17] the authors show that a coherent 
(non-monotonic) path traversal that uses Grover's fixed 
point method for this case can be implemented with an 
improved cost C(Llog 2 L/VT). 



V. RELATION TO OTHER WORK 

It has been noted previously [TO] [TJ EH HZ] that 
the projective-measurement operations Mi can be simu- 
lated using Kitaev's phase estimation algorithm [TT] in 
the discrete-time case. This requires implementing uni- 
taries Ui — e~ lH ^ controlled on r ancillary qubits initial- 
ized in the equal superposition state. Then the inverse of 
the quantum Fourier transform is applied to the ancillary 
qubits, and a projective measurement on the computa- 
tional basis of the ancillae is performed [see Fig. [It a)]. 
The phase estimation algorithm needs to resolve the de- 
sired eigenphase from other eigenphases to be able to 
project the state of the system into the desired eigenstate. 
This requires 2 r 6 il(l/A(l)) uses of controlled- Ui's for 
constant error. The error per step has to be small. If one 
of the high-confidence versions of the phase estimation 
algorithm |28] is used, the overhead to achieve error e is 
logarithmic in 1/e. The overall cost is then similar to 
that of the randomization method when T is restricted 
to be positive. 

Interestingly, the phase-estimation-based algorithm 
produces the same effect on the system as the random- 
ization method if we sample the evolution time from the 
uniform distribution on an interval. This is because the 
phase estimation ancillary qubits can be traced out af- 
ter each step. As a result, the inverse quantum Fourier 
transform can also be dropped. Consequently, the co- 
herence in the state of the ancillary qubits, initialized 
in the equal superposition state, plays no role and these 
qubits can be replaced by classical bits, each being or 1 
with probability 1/2. This equivalence was also studied 
in Ref. [2J]. We illustrate it in Fig. [TJb). 
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FIG. 1: (a) Phase estimation algorithm. At the end of the 
algorithm, the top r-qubit register encodes a r-bit approxima- 
tion to an eigenphase of Ui on readout. It is initialized with 
Hadamard gates to an equal superposition state. A sequence 
of 2 r — 1 controlled Ul operations is applied, and the first reg- 
ister is measured after an inverse quantum Fourier transform. 
If the measurement outcome approximates an eigenphase of 
Ui, the second register (system) is approximately projected 
onto the corresponding eigenstate. (b) Randomized evolu- 
tion. If the phase estimation algorithm outcome is ignored, 
the overall effect is equivalent to the one induced by initial- 
izing a set of r bits (first register) in a random state j, with 
j £ [0..2 ? — 1], and by acting with Uf. Double lines indicate 
classical information. 



Repeating the phase estimation algorithm n times is 
equivalent to randomizing with the sum of n independent 
uniform distributions. This was considered in Example [3] 
Sec. |II C[ The unwanted coherences reduce exponentially 
in n. 

There are previously noted relationships between the 
Zeno effect and coherent evolutions similar to the con- 
tinuous or discrete evolutions used in the randomization 
method. For example, the effect of a strong interaction 
with another system, such as might occur in the coupling 
to a measurement apparatus, is to restrict the natural 
Hamiltonian to the eigenspaces of the interaction [29] , 
The suppression of coherent transitions by randomiza- 
tion with the interaction Hamiltonian would have a sim- 
ilar effect. A discrete version of this observation relevant 
to the analysis of dynamical decoupling was considered 
in 0U]. 

There is a relationship between the way in which inter- 
actions are averaged away in dynamical decoupling, par- 
ticularly randomized dynamical decoupling |311 132] , and 
how transitions between the adiabatic path and the other 
eigenstates are suppressed in the randomization method. 
The relationship can be made explicit by changing to an 
s-dependent frame in which the Hamiltonians H(s) are 
diagonal. In this frame, the transitions show up explicitly 
due to the frame changes with s. Strong or randomized 
evolution under H(s) suppresses these transitions by av- 
eraging them to zero. Dynamical decoupling typically 
uses operators that have stronger averaging effects. 

A feature of the randomization method is the use of 
phase decoherence to ensure a more efficient transfer 
to a state of physical or computational interest. There 
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are other ways in which decoherence can play a role in 
preparing states for quantum computing. Early examples 
proposed the synthesis of decoherence free subspaces from 
an environment-induced quantum Zeno effect [33 , 3Hl35|. 
The use of decoherence to decrease the mixing time of 
quantum walks was proposed in Ref. |36j . A related 
phenomenon has been studied in the context of energy 
transfer [371 [351 IMJ SO] as realized in certain biological 
molecules. More generally, it may be that decoherence 
or thermal noise can enhance the success probability in 
adiabatic quantum computing [UJ H21 23J [H] . Note that 
the required thermal noise is different from the phase- 
decoherence associated with the randomization method 
in that it has the potentially desirable effect of trans- 
ferring population to lower-energy eigenstates of the cur- 
rently active Hamiltonian. Whether the requirements for 
effective exploitation of this situation can be met in re- 
alistic devices is not clear. 

Finally, the engineering of a dissipative process using 
feedback techniques to stabilize a desired quantum state 
has been studied extensively in the quantum control lit- 
erature j45j|46j|47l|48]. Reference [49] advances that the 
final state of AQC can be made the unique steady state 
of the dissipation, even without feedback. In essence, 
the process' Lindblad operators encode the gates of the 
quantum computation and appropriate updates to a log- 
ical clock register. Again, the necessary dissipation re- 
quires more than the phase decoherence realized by the 
randomization method. 



VI. CONCLUSIONS 

We have described a method for state preparation in 
the spirit of AQC, but based exclusively on randomized 
evolutions. The idea is to perform a discrete sequence of 
projective measurement operations onto the desired (in- 
stantaneous) eigenstate of a given Hamiltonian or uni- 
tary path. These operations are induced via evolution 
randomization, which realizes the necessary decoherence 
in the eigenbasis. We bound the residual coherences after 
the randomization in terms of the characteristic function 
of the random time. 

We obtained the following exact bounds on the dephas- 
ing achieved by randomized evolutions: First, to induce 
enough decoherence, the average evolution time per step 
scales with the inverse of the minimum absolute value of 
the spectral (or eigenphase) gap. Second, repetition of 
the randomization reduces the coherences exponentially 
in the amount of repetitions. Third, if negative-time evo- 
lutions are implementable with constant overhead, loga- 
rithmic factors depending on the error can be reduced to 
constant factors, even for discretized evolutions. Fourth, 
for non-negative evolutions and if only a lower bound on 
the absolute value of the gaps is known, the logarithmic 
overhead is unavoidable. 

We show that the complexity of path traversal algo- 
rithms is best expressed in terms of the path length L. 



The explicit dependence of the complexity on L can be 
very helpful when L does not depend on the gap. This 
happens, for example, in the Hamiltonian version of an 
algorithm for unstructured search, where we showed that 
a simple choice of step size and random time distributions 
rotates into the solution state with probability 1/2. One 
further advantage of the path-length formulation is that 
we do not require the relatively strong differentiability 
requirements on H as in the proofs of the adiabatic con- 
dition [13] with explicit bounds as needed for AQC. 

Another case where L does not depend on the gap is 
in the quantum simulated annealing algorithm, which we 
also analyzed. This algorithm provides a quadratic quan- 
tum speed-up in terms of the gap with respect to classi- 
cal simulated annealing implemented via Markov Chain 
Monte Carlo methods. The path is determined by an 
annealing schedule in which a parameter (3, related to 
the inverse temperature of a classical system, is slowly 
increased in equal-size steps. The quantum simulated 
annealing algorithm allows us to reach the optimal con- 
figuration in time 0(1/VT) for constant probability of 
success and path length, with T being the minimum gap 
of the stochastic matrix (and the corresponding Hamil- 
tonian) along the path. The improved randomization 
methods given here remove a logarithmic factor for the 
version of the algorithm given in Ref. [12 . 

The similarities of the randomization method with 
AQC are clear: A typical instantiation (choice of evo- 
lution times) of the randomization method is, with high 
probability, an approximation to an adiabatic path. We 
find, as is often the case, that it is easier to prove er- 
ror bounds for random instances than for the worst case. 
Whether the existence statement can be "derandomized" 
efficiently is still an interesting question. 
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Path length. For states \<fii) and \4>2), let 
8(|</>i) , |02)) = arccos(|(02|0i)|) be the angular distance 
between the states. We assume that the \ip(s)) form a 
projectively continuous path, s G [0,1]. The length is 
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given by 



l = sup^e(|^(sfc+i)) , \i> Sh )), 

(s k ) k 



(21) 



where the ordered sequences (s^) subdivide [0,1]. Note 
that the expression in the limit depends monotonically 
on the (sfe), increasing in the refinement order. If \ip(s)) 
is differentiable, the expression in Eq. (21 1 reduces to the 
one m Eq. ^. 

Uniform parametrization. Let L(l) be the length 
of the path \ip(l')) for < I' < I, defined as in 
Eq. (21 1. Suppose that L(l) is Lipschitz continuous so 

that w(h,h) = ^P h <i'<i"<i 2 (L( l ") ~ UiW" ~ 
is finite. Note that if L is differentiable, one can take 

w(h,h) = su Ph<i<i 2 ^sP~ ■ In particular, if |0(/)) is 
differentiable, u(l\,l2) = sup /l<;< ; 2 ||9i|V>(Z))|| works [see 
Eq. P 



We obtain: 



Lemma 7. The squared overlap \ (I + 5)1^(1)) \ 2 can be 
bounded by 



2x2 



\W + S)m))\ >l-u(l,l + S) J 6 



Proof. We have 

\w+s)\m)\ 2 



= cos(e(W + 5)} Am))) 2 

> i-e(m+s)),\m)) 2 

> 1 -(£(/ + S) -L(l)) 2 

> 1-uj(1 + 5,1) 2 5 2 (22) 



To take advantage of Lemma[T] it helps to parametrize 
the path with an s(l) for which Lo^h,^) is as uniform as 
possible. For this purpose, define s(l) = inf{s : L(s) > 1} 
for < I < L, where the length L(s) is the length of 
the path \tp(s')), < s' < s. The function s(l) is not 
necessarily continuous. 

Continuity of states and finiteness of L implies con- 
tinuity of L(s). This can be shown as follows: Sup- 
pose that L(s) is not continuous at s. Then ei- 
ther sup 5>0 L(s — S) < L(s) or inf^x) L(s + S) > 
L(s). Consider the first case. We have L(s) — 
limsup, >0 (6(|V(s)) , |^(s - 5))) + L(s - 5)). The in- 
equality implies that lim sup (5>0 Q(\ip(s)) , 1-0(5 — <>))) > 
0, contradicting continuity of \tp(s)). For the second 
case, s < 1. Define L(si,S2) as the length of the 
path from \tp(si)) to \ip(s2))- It can be seen from the 
definition, monotonicity in the refinement order of the 
term in the limit of the definition, and from projec- 
tive continuity of \ip(s)) that L(s, 1) = sup (5>0 L(s + 
8,1) + 0(|^(a)),|V(* + O» = sup 5>0 L(s + ,5,1) and 
L(s,l) = L(s + 6,1) + L(s,s + S). It follows that 
inf^x) L(s, s + S) =0. The observation now follows from 
L(s + S) = L(s) + L(s, s + 5). 

We define L(l) as the length of the path \tp(l')) = 
|V>(s(Z'))) for </'</. We show that 1(1) = L(s(l)) = I. 



The second inequality follows from continuity of L and 
the definitions. From the definition of path length and 
since any subdivision (l k ) of [0, Z] corresponds to a sub- 
division (s(Zfc)) of [0,s(/)L L(l) < L(s(l)). To show 
the reverse inequality, let s = s(L(s)). Then s < s 
and Q(\ip(s)} ,\i>(s))) = 0. Hence for all s' G [s,s], 
\ip(s')) oc \ip(s)) (that is, the two states are projectively 
identical). Consequently, the right-hand side of Eq. (21 1 
is unchanged if we replace the Sk by Sk ■ Since the Sk are 
in the range of / i— > s(l), we can choose lk = i(sfc) to 
show that the defining supreme for L(l) and for L(s(l)) 
are the same. 

By the previous paragraph, lo(Ii,12) — 1 for the 
parametrization s(l). We therefore refer to s(l) as the 
uniform parametrization. 

Proof of Theorem [l] Let fi be the probability distri- 
bution of T. For any £ , 

(nf-Mf)(\m)m)\) 
= njmi))mm 

- ik ^\\m)(ui)\v mt d^t) 



(23) 



e^dn{t)\m){i>M 
$(^)ivi(0)(^(0l • 



We assume without loss of generality that p is pure, p = 
\4>){<j>\. Write 



ei|#)) + 5>|^(Z)) 



(24) 



tion TZf leaves S invariant, and we can choose £ = TZf 



Let S be the subspace orthogonal to |i/>(Z)). The opera- 
tion TZf leaves S invarif 
in that subspace. Then 

WiMf-nfm (0|)|| tr 

Kf (^c lC *|0(O)(^(OI+h-c 
\j">i 



J>1 



This is the trace norm of a matrix having 



(25) 



± /> m^)c lC ;\ 2 




(26) 



as the only non-zero eigenvalues. Because of the normal- 
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ization, \cx\ 2 Y, 3 >x \c 3 \ 2 < 1/4. Thus 

\\(Mf~nfm (0|)|| tr 

= 2 Stic; 
Vi> 1 



* 12 



<sup|$( Wj -)|2 4 /£><$| s 

V j>i 

<sup|$(^)| . 



(27) 



Proof of Lemma [3] We start with any smooth even 
function h of compact support in (—1/2,1/2). This 
implies that its inverse Fourier transform h is real and 
all its moments are bounded since 



\(X»)\ 



+oo 



h(x)x n dx 





d n h(0) 




duj n 



< oo 



(28) 



We define the characteristic function <!>! to be propor- 
tional to the convolution of h with itself, 



oc (h * h)(uj) . 



(29) 



We normalize such that $1(0) = 1. By construction, 
the inverse Fourier transform of $1, denoted by fx, is 
positive, normalized to 1, and rapidly decaying, as de- 
sired. To accommodate arbitrary spectral gaps A > 0, we 
rescale the characteristic function as <&a(w) = 3>(w/A), 
which has support in (—A, A). Its inverse Fourier trans- 
form is a probability density function /a(^) = A/i(Ai). 
The cost of randomization with /a is 

/+00 
\t\f A (t)dt 
-OO 

/-j-00 
\t\h(At)dt 
-OO 



|T|>i 



(30) 



where (|T|)i is the cost of randomization with fx, and 
is independent of A. It follows that (|T|) A £ 6(1/A), 
which is optimal. 

Proof of Lemma [4] Consider a probability density 
/ with characteristic function $ of support in (—A, A), 
where A < ir. Consider J^k^i^ + 2lrk ) = * C)(oj) 
where C(oj) — J2k ~ 2nk) is a comb. As a dis- 
tribution, C{uj) is the Fourier transform of the comb 
C{t) = £ fe <5(i ~ fc VO")- See > for example, Sec. 2.4 
of [50] ■ Using the rules for convolution under the in- 
verse Fourier transform, we find that the distribution 
(f-C)(t) = Efc/( fc )<50- fc ) has Fourier transform $* C. 
Because ($ * C)(0) = $(0) = 1, it follows that f(k) is a 
probability distribution with the stated properties. 

Proof of Theorem [3j For (T) infinite, there is noth- 
ing to prove. So assume (T) is finite, which implies that 



the characteristic function is differentiable. Suppose first 
that T has a square-integrable probability density f(t). 
The characteristic function is then a "Hardy function" of 
class H 2+ as defined in Ref. [ST], pg. 162. By noting that 
for a > 0, u> 1— > Q(auj) is also Hardy, the proof of Thm. 2 
on pg. 166 of Ref. [51] shows that 



+ °° loff|$(cnOI 



OO 



I + 7 2 



(31) 



where $ has been analytically extended to the upper 
half plane. The analytical extension of $ is obtained 
by using complex lu in the Fourier transform. Conse- 
quently, d$>(z)/dz is the Fourier transform of 1 1— > itf[t), 
where defined. In particular, |d<i>(z)/(iz| is bounded by 
(T) for z = i/3 with /3 > 0. Since $(0) = 1, we have 
log |$(m)| > log(l — a(T)). The integral of the inequal- 
ity in Eq. ( |31[ ) can be related to the desired supremum 
as follows: 

'+°° log|*(a7)| 



I + 7 2 



-d~j 



< a 



+o ° log \m l 

■ 00 " + 7 

+A 



A 



log |*(7)l . 

q,Z _|_ 



+ alog( sup I $(7) |) 
M>A 

< log( sup |*(7)|)(ir 
ItI>A 



1 



| 7 |>A a +7' 

2 arctan(A/a)) 



:d7 



(32) 



To drop the first summand in the last step we used the 
fact that |3>(7)| < 1 because <& is the characteristic func- 
tion of a probability distribution. We now let a — > + 
and combine with the earlier inequality to get, to first 
order in a, 

- na{T) < 21og( sup |$(7)|)a/A , (33) 
M>a 

which gives e~ A ^ T ^ < sup| 7 | >A |3>(7)|- 

Now consider arbitrary positive T with (T) < 00, 
and with probability distribution /1. Let Sg be uni- 
formly distributed between and S. The probabil- 
ity distribution of T + S has cumulative distribution 
F(x) — Jq min(l,(x — y)/S)dfi(y), which is differen- 
tiable. The corresponding probability density is given by 
fj,([x — S, x])/S = f g dfj,(y)/S and is square integrable 
because 



K[y - 5 ,v\) 2 dy < 



fj>([y- s > y\) d v 

V 



y-8 
z+S 



dfi(z)dy 



dydfi(z) 
5d[i(z) = S . 



(34) 
(35) 
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Thus T+ S is subject to the bound of the Theorem. The 
characteristic function of T + S is given by $>(u))ss(u>), 
where Ss(u>) is the characteristic function of S$. The 
function ss{u>) converges uniformly to 1 on bounded in- 
tervals as 6 — > + . It follows that the desired bound 
applies to arbitrary positive T. 

Proof of Lemma |5} Without loss of generality, the 
phases of \ip(s)) are geometric. Because A(s) > and 
H(s) is differentiable, it follows that \ip(s)) is differen- 
tiable. From the eigenvalue equation 



we get 



H (s) | V(s)> = E( 3 )Ms)) , 



d s H(s)\ij(s)) + H(s)\d s ij(s)) = 
d s E( S )\yj( s )) + E(s)\dMs)) . 



(36) 



Denote by \ipj{s)), j € {2, 



(37) 

d}, the j-th eigenstate 



of H(s), orthogonal to \ip(s)), and with eigenvalue Ej(s). 
We obtain 

Because the path \ip(s)) is geometric, (■0(s)|9 s i/'(s)) = 
for all s. This gives 

iiw.)>i--E "^^y 



J>2 



J>2 



1 

A (7) 2 



<X7^(^)|(^( S )) 2 |^( S ))< 



d s g(s) || 
A(s) 2 
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